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I. INTRODUCTION 

The paper is devoted to the problem of constructing a special class of symmetries for boundary value prob¬ 
lems (BVPs) in mathematical physics, namely renormalization group symmetries (hereafter referred to as RG- 
symmetries). 

Symmetries of this type appeared about forty years ago in the context of the renormalization group (RG) con¬ 
cept. This concept originally arose jl] [|] in the ’’depth” of quantum field theory (QFT) and was connected with a 
complicated procedure of renormalization, that is ’removing of ultra-violet infinities’. In QFT renormalization group 
was based upon finite Dyson transformations and appeared as a continuous group in a usual mathematical sense. 
It was successfully used for improving approximate perturbation solution to restore a correct structure of a solution 
singularity. 

In the seventies, it was found that the RG concept was fruitful in some other fields of microscopic physics: phase 
transitions in large statistical systems, polymers, turbulence, and so on. However, in some cases, following Wilson’s 
approach || to spin lattice, the original exact symmetry underlying the renormalization group notion in QFT was 
changed to an approximate one with the corresponding transformations forming a semi-group (not a group as in 
the QFT case). Here, in this paper, by RG-symmetry we mean the original exact property of a solution - as 
it was formulated in Refs. ||,Q] (see also J(|) by N. Bogoliubov and one of the present authors. Thus, by RG- 
symmetry we mean a symmetry that characterizes a solution of a BVP and corresponds to transformations involving 
both ’’dynamical” (i.e., equation) variables and parameters entering into a solution via equations and boundary 
conditions. 

For a simple illustration we consider some BVP that produces a family of solutions. The simplest variant of RG 
transformation is given by a simultaneous one-parameter point transformation 

T a : {x -> x' = x/a , g -> g' = G(a,g) } , G(l,g) = g (1) 

of a dimensionless ’’coordinate” x and a one-argument ’’characteristic” (e.g., initial value) g of each solution, the 
quantity of a direct physical interest. The transformation function G(x,g), which depends upon two arguments [[zj, 
should satisfy the functional equation 

G(x, g) = G{x/a , G(a , g)), (2) 

that guarantees the group property T a ■ Tj, = T a b fulfillment. 
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The functional equation (|j) and transformation (|l]) arise, for example, in the massless QFT with one coupling. In 
that case x = Q 2 / p 2 is the ratio of a 4-momentum Q squared to a ’’reference momentum” p squared, g is the coupling 
constant and G is the so-called effective coupling. 

Later on the -type symmetry underlying the renormgroup invariance was also found in a number of problems 

of macroscopic physics like, e.g., mechanics, transfer theory, hydrodynamics and a close relation of RG-symmetry to 
the notion of self-similarity was established | £ |uj . 

The infinitesimal form of transformation (1) can be written down as a differential equation 


RG = 0, with R = xd x — (3(g)d g , (3(g) = 


dG(a\g) 


da 


(3) 


where R is the infinitesimal operator of RG-symmetry (or, simply, RG-operator) with the coordinate (3(g) defined by 
the derivative of the function G. 

Therefore, instead of relations (® and ([|), RG transformation can also be introduced by means of an RG-operator. 
And vice versa, being given an RG-operator one can reconstruct the functional equation for the solution G with 
the help of the characteristic equation for (j||). Moreover, for a given (3 -function or, in other words, for the given 
RG-symmetry, one can get an explicit expression for the invariant of the group transformation G(x, g) by solving the 
corresponding Lie equations |ll[] 


dx' 

x' 


dg' 

W) 


da 

a 


(4) 


with the initial conditions x'\ a —i = x, g '| a =i = g. 

Along with (^]) a different form of the invariance condition for the function G(x, g) is often employed 


— G(x/a,G(a,g)) 


= 0 . 


(5) 


a— 1 


Equation (^|), reflecting the invariance of G under the RG transformation can be treated as a vanishing condition for 
the coordinate ae of the RG operator (|j) in the canonical form (l^] 


R = eedc , ae = xG x — (3(g)G g = 0, 


( 6 ) 


identically valid on a particular BVP solution G = G(x,g). 
At the same time, the relation 


RS(x,g) = (xd x - (3(g)d g ) S(x,g) = i(g)S(x,g) (7) 

corresponds to the function S(x,g) that is a covariant [^3| of the RG transformation. In QFT case, this relates, e.g., 
to a propagator amplitude (see, Refs. Jfj] and H). Here, 7 (g) is known as the anomalous dimension of S(x,g). 
Generally, the differential equation akin to (f|) 


xfx ~ (3(g) f g = 0 . ( 8 ) 

states an invariance of a function / under the RG transformation (Jl). Its solution f(x,g) = F(G(x,g)) precisely 
corresponds to the same property emphasized by (|j). 

In a particular case, when the function G is linear in the last argument, G ~ kg, equation (||) defines a solution 
that has a power x dependence, i.e., G(x,g) = gx k with k being an arbitrary number. Then, equation (^|), takes a 
form of power scaling (or power self-similarity) transformation 

x' = x/a , g' = ga k , 

that is well-known in mathematical physics and widely used in the problems of hydrodynamics of liquids and gases. 

Therefore, transformation (j|) can be considered |I(J as a functional generalization gx k —> G(x,g) of a usual (i.e., 
power) self-similarity transformation. One can refer to it as to functional self-similarity transformation: this term 
was first introduced in H as a synonym of the RG transformation as defined above. 

It is widely known that in QFT, as well as in other mentioned fields of theoretical physics, the RG method allows 
one to improve the perturbation theory results and to simplify the analysis of a singular behavior of a solution which 
becomes scale-invariant in the vicinity of a singularity. The latter reminds a situation, which is typical of asymptotic 
analysis of solutions of differential equations (DEs): long-time asymptotics demonstrate self-similar regimes [[Td||. 
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Hence, it looks natural to use the RG methods to study strong nonlinear regimes and to investigate asymptotic 
behavior of physical systems described by DEs. We have no possibility to discuss this in detail here and would limit 
ourselves to mentioning some successful attempts of using the RG ideas in mathematical physics. 

To our knowledge, the very first results in this field were obtained about a decade ago by two of the co-authors (k| 
of this paper by applying RG ideas to a problem of generating of higher harmonics in plasma. This problem, after 
some simplification, was reduced to a couple of partial DEs with the boundary parameter - solution ” characteristic” - 
explicitly included. It was proved that these DEs admit an exact symmetry group similar to that, defined by Eq.(^|). 
The group obtained was then utilized to construct the desired nonlinear solution of the BVP. This approach has 
further been developed (see jT(| and references therein) and we shall describe it in some detail in the next Section. 

The methods of QFT RG were exploited by Goldenfeld, Martin and Oono with co-authors (Urbana group) ( see, 
e.g., Refs, [jnj and [0) to find asymptotics of the solutions of parabolic-type nonlinear differential equations, that 
describe a variety of physical phenomena, such as groundwater flow under gravity, shock waves dynamics, radiative 
heat transfer and so on. As an auxiliary tool, they used the concept of intermediate asymptotics first introduced by 
Barenblatt and Zeldovich 0 - see also the review monographs 0 and [0 . In this way, the Urbana team was able 
to determine values of exponents in the ratios of invariants forming arguments of self-similar solutions. 

Later on, with the goal to a global asymptotic analysis they developed and illustrated, by numerous examples, the 
”perturbative renormalization group theory” (see |pl|| and references therein) that exploited the form of the invariance 
condition akin to that (|j) used in QFT. The geometrical formulation of the perturbative RG theory for global analysis 
was presented by Kunihiro 0) on the basis of a classical theory of envelopes. 

On the other hand, Bricmont and Kupiainen |2f| j2f| attracted RG ideas for nonlinear DEs analyzing in a bit 
different manner. They used an iterative set of rescalings borrowed from the Wilson version of renormalization group, 
that is semi-group. On basis of that RG-mapping procedure they succeeded in proving the global existence and 
detailed long time asymptotics for classes of nonlinear parabolic equations. 

Generally, the procedure of revealing RG transformations, or some group features, similar to RG regularities, in any 
partial case (QFT, spin lattice, polymers, turbulence and so on) up to now is not a regular one. In practice, it needs 
some imagination and atypical manipulation (see discussion in jfi . 001 ) ’’invented” for every particular case. For 
example, the above described RG methods applied to asymptotic analysis of differential equations were based on the 
a priori assumption of the existence of some scaling transformations or on the invariance condition of an approximate 
solution. By this reason, the possibility to find a regular approach to constructing RG-symmetries is of principal 
interest. In this paper we give an account of our efforts for creating a possible scheme of this kind in application to 
physical systems that are described by DEs. The leading idea in this case is based on the fact that symmetries of 
such systems can be found in a regular manner by using the well-developed methods of modern group analysis. 

The paper is organized as follows: in Section II we describe the general scheme of constructing RG-symmetry for 
a BVP. It appears, that the implementation of this scheme strongly depends on the mathematical model used and 
on the form of boundary conditions. As a result, different approaches to finding RG-symmetries are possible, and 
these are illustrated by examples in the following five sections. In Section III RG-symmetries are calculated using the 
classical Lie symmetries. In Section IV RG-symmetries are obtained on the basis of Lie-Backlund symmetries. In the 
next two sections RG-symmetries are found when boundary conditions are presented either in the form of a differential 
constraint (Section V), or in the form of an embedding equation (Section VI). In Section VII one more approach to 
RG-symmetries constructing is presented which is based on an approximate group symmetry. In conclusion, we make 
a summary of the approach and discuss some further applications. 


II. APPROACH TO CONSTRUCTING RG-SYMMETRIES 

First of all, we emphasize that the desired regular approach to constructing RG-symmetries turns out to be possible 
for those mathematical models of physical systems that are based on differential or, in some particular cases, integro- 
differential equations. The key idea uses the fact [0,0] that such models can be analyzed by algorithms of modern 
group analysis. 

The proposed scheme comprises a sequence of the four steps. 

I. A specific manifold (differential, integro-differential, etc.) should be primarily constructed. This manifold that 
will be referred to as renormgroup manifold (RG-manifold) generally differs (see below) from the manifold given by 
the original system of DEs. 

II. The second step consists in calculating the most general symmetry group Q admitted by the RG-manifold. 

III. The restriction of the group Q on the desired BVP solution (exact or approximate) constitutes the next 
step. The group of transformations thus obtained (renormgroup) is characterized by a set of infinitesimal operators 
(RG-operators), each containing the solution of a BVP in its invariant manifold. 
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IV. The last, fourth step implies utilization of RG-operators to find analytical expressions for solutions of the BVP. 

Being formulated in a concise form these steps deserve further comments. 

Comment to I. In the scheme described above the first step, namely constructing the RG-manifold, is of fundamental 
importance. The form of its realization depends both on a mathematical model and on a form of a boundary condition. 

Here we show the following different approaches to RG-manifold constructing: 

la. In the first, more simple case the RG-manifold, as usual in classical group analysis, is presented by a system of 
basic DEs with the only substantial difference: parameters, entering into a solution via the equation and boundary 
conditions, are included in the list of independent variables. 

l b. Another approach to constructing the RG-manifold implies an extension of a space of variables involved in 
group transformations, for example, by including differential variables of higher order and nonlocal variables. It means 
that in this case Lie-Backlund transformation groups and nonlocal transformation groups should be invoked mi¬ 
le. In the third approach the procedure of construction of RG-manifold is based on the invariant embedding method 

po| . Here RG-manifold is given by a system of equations that consists of original DE and/or embedding equations 
which correspond to the BVP under consideration. 

l d. The fourth approach to some extent is similar to the previous one. In this event boundary conditions are 
reformulated in terms of a differential constraint which is then combined with original equations to form the desired 
RG-manifold. 

l e. The last approach utilizes approximate transformation groups. Here, the RG-manifold is given by a system of 
DEs with small parameters and can be analyzed by perturbation methods [[il| . 

Comment to II. Searching the symmetry of RG-manifold is the main problem of the second step. The term 
’’symmetry” as used in the classical group analysis means the property of a system of DEs to admit a Lie group 
of point transformations in the basic space of all independent and dependent (differential) variables entering these 
DEs. The Lie calculational algorithm of finding such symmetries is reduced to constructing tangent vector fields with 
coordinates, that are functions of these basic group variables and can be defined from the solution of an overdetermined 
system of DEs, named as determining equations. In modern group analysis different modification of a classical Lie 
scheme are in use (see, e.g. |2ti| . [i2f and references therein). If the problem of finding symmetries for a given system of 
DEs (RG-manifold) is solved, then the result is presented in the form of Lie algebra of infinitesimal operators (also 
known as group generators), which correspond to the admitted vector field. In what follows these operators will be 
denoted by X. 

Comment to III. The goal of a group restriction is the construction of a transformation group with a tangent vector 
field (point, Lie-Backlund, etc.) infinitesimal operators of which (hereinafter referred as R) contain the desired BVP 
solution in an invariant manifold. This means, that the coordinate of the canonical operator of RG-symmetry vanish 
on the BVP solution and on its differential consequences. 

Mathematically, the procedure of a group restriction appears as a ’’combining” of different coordinates of group 
generators X admitted by the RG-manifold. The vanishing condition for this combination on a solution of the BVP 
leads to algebraic equalities that couple different coordinates and give rise to desired RG-symmetries. In a particular 
case, when RG is constructed from a Lie group admitted by the original system of DEs, it turns out to be a subgroup 
of this group and a solution of the BVP appears as an invariant solution with respect to the point RG obtained 
(compare with |ll]]). In the general case, not only Lie point group, but Lie-Backlund groups, approximate groups, 
nonlocal transformation groups, etc. (see, e.g. [^9)), are also employed as basic groups which are then to be restricted 
on the solution of a BVP. 

Comment to IV. A technique for constructing group invariant solutions corresponding to a symmetry group when 
its infinitesimal operators are known has been detailed in various monographs (see, e.g., ). Therefore, the 

final step is performed in a usual way and needs no specific comments. 


Before proceeding any further, we make a short review of results, that were obtained on the basis of the formulated 
scheme. The first application of RG-approach to a particular problem of laser plasma was announced in | jl5| . This 
problem, namely the problem of a nonlinear interaction of a powerful laser radiation with inhomogeneous plasma, 
has been detailed in subsequent publications [fi7i|.fS3 341. A mathematical model was given by a system of nonlinear 


DEs for components of electron velocity, electron density and the electric and magnetic fields. The presence of small 
parameters (such as weak inhomogeneity of the ion density, low electron thermal pressure and small angles of incidence 
of a laser beam on plasma surface) in the initial system of equations provided a way to constructing RG-manifold 
using (Ie) approach, based on approximate group methods. The desired RG-symmetry appears as Lie point symmetry 
that takes account of transformations of a boundary parameter (common to (la) approach), which is related to the 
amplitude of the magnetic field at a critical density point. RG-symmetry obtained made it possible to get the exact 
solution of original equations, that was then used to evaluate the efficiency of harmonics generation in cold and hot 
plasma (see [pT|). 
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The advantageous use of the RG-approach in solving the above particular problem gave promise that it may work 
in other cases. This was illustrated in |35|| by a series of examples for different BVPs. Various methods of constructing 
RG-symmetries were described, based on the use of point symmetries, approximate symmetries, embedding equations 
and transformations of Fourier components. As is shown in |j^], different formulations of BVPs give rise to various 
methods of finding RG-symmetries. Thus, further development of the scheme was concentrated on analyzing these 
methods. 

The first one was concerned with the initial value problem for the modified Burgers equation with parameters of 
nonlinearity and dissipation included explicitly. This example f36| yielded a detailed illustration of the method of 
constructing RG-symmetries when a basic RG-manifold is given by an original DE with parameters included in the 
list of independent variables (la approach). It was argued that the exact solution can be reconstructed from the 
perturbative solution with the help of any of the admitted RG-symmetry operators which form an eight-dimensional 
algebra. Two illustrative examples were given, dealing with perturbation theory in time and in nonlinearity parameter. 

To demonstrate the method of constructing Lie-Backlund RG-symmetries that uses (lb) approach, the initial 
value problem for a linear parabolic equation was considered in ||27| ]. It was shown that appending Lie-Backlund 
RG-symmetries to point RG-symmetries extends the algebra of RG-symmetries up to an arbitrary order. 

The same mathematical model was also employed within the (Ic) approach when the boundary condition is described 
by a differential constraint |37| . It was found that RG-symmetries obtained can not be reduced to point RG-symmetries 
which arise from the (la) case. However, some of them can be reformulated in terms of RG-symmetries previously 
found in (la) approach while the others can be constructed from the Lie-Backlund symmetries of basic equations in 
view of the given differential constraint. 

An idea of the (Id) approach to constructing RG-manifold based on the invariant embedding method was realized 
in p5| ] for ordinary DEs. Here, embedding equations can be treated as a specific form of a differential constraint, 
that takes boundary data into account. The method of finding RG-symmetry using the (Id) approach proves to be of 
particular interest for the first order ordinary DEs when using (la) approach faces standard problems in calculating 
point symmetries admitted by RG-manifold. Provided the embedding equation has the form of a first order DE of 
an evolutionary type there appear no difficulties in group analysis of a joint system of the basic and the embedding 
equations. 

Worthy of mention is an example that demonstrates the utilization of RG-symmetries to constructing solutions of 
the BVP for a system of two first-order partial DEs that describes the propagation of a laser beam in a nonlinear 
focusing medium (3^| (|i| . It was revealed that RG-symmetries are related to formal symmetries that are constructed 
in the form of infinite series in medium nonlinearity parameter. For a specific form of boundary data infinite series 
are truncated with RG-symmetries presented by finite sums. Generally, for arbitrary boundary data this is not the 
case and in that event a finite sum describes approximate RG-symmetry for small nonlinearity parameter. Based on 
(la), (lb) and (Ie) approaches both point and Lie-Backlund (exact and approximate) RG-symmetries were obtained 
and then used to find an analytical solution of the problem. 

To clarify the idea of constructing RG-symmetries several examples are given below which demonstrate different 
approaches to the problem. To gain better understanding of these approaches, a simple mathematical model is used 
[tfjf . This model corresponds to BVP for a system of two first-order partial DEs that were studied by Chaplygin 
in gas dynamics 


v t + vv x - aip(n)n x = 0 , n t + vn x + nv x = 0 ; 
n(0, x) = V ( x ), n(0, x) = N(x ), 


(9) 


where <-p{n) is an arbitrary function of n and a is a nonlinearity parameter. Despite its simplicity, this mathematical 
model has a wide field of application and was used to describe various physical phenomena (in the so-called quasi- 
gaseous media pdl ). In such a case the physical meaning of variables f, x, v and n may differ from that in gas 
dynamics. For example, in nonlinear geometrical optics, t and x are coordinates, respectively, along and transverse to 
the direction of laser beam propagation, v is the derivative of eikonal with respect to x, and n is a laser beam intensity. 
In this case, functions V and N characterize the curvature of the wave front and the beam intensity distribution upon 
the coordinate x and the entrance of a medium t = 0 . 

Along with ([jj) another form of basic equations will be used 

ar v - ( n/ip(n))xn = 0 , \v + x n = 0 . (10) 

These linear equations for new variables r = nt and \ = x — vt results from (|(jj) under hodograph transformations. 
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III. RG AS LIE POINT SUBGROUP 

This section presents an illustration of the method of constructing RG-symmetries when a basic RG-manifold is 
given by the original DEs with parameters included in the list of independent variables. Boundary conditions are taken 
into account while restricting the group admitted by RG-manifolds up to the desired RG on the exact or approximate 
solution of a BVP which thus appears as an invariant solution with respect to any of RG operators obtained. 

3.1. First we shall consider a particular case of equations (j9) when the nonlinearity parameter a is equal to zero: in 
application to optical equations discussed above this means that nonlinear effects are neglected. Then by introducing 
a new variable v = su, the system of equations (|)j) is rewritten in the following form 

Ut + euu x = 0 , nt + eun x + enu x = 0 ; ( 11 ) 

u(0,x) = U(x), n(0,x)=N(x). ( 12 ) 

The continuous point Lie group admitted by the differential manifold © (RG-manifold) is given by the infinitesimal 
operator (a general element of Lie algebra) with six independent terms 


6 

X = ^dt + £ 2 d x + £ 3 d e + ifd u + ri 2 d n = ^ Xj , (13) 

i =1 

Xi = (1/e) AJ 1 dt + (J 1 + u A J 1 ) d x - nJld n , A J k = ( etJ* - J k ) , 

X 2 = (1 /n) J 2 {d t + eud x ) , X 3 = nJ 3 d„ , X 4 = J 4 {-td t + nd n + e<9 e ) , 

X 5 = AJ 5 D + J 5 {£td x + d u ) , X e = — (1/n) J 6 D , D = (td t + £utd x — nd n ) . 

Coordinates £ and rj of this infinite-dimensional group operator depend upon five functions J l {x-, u, £), * = 1 , 2 ,3, 5, 6 
which appear as arbitrary functions of their arguments x = x ~ vt, u and e. The sixth one, J 4 , that enters the operator 
which describes group transformation of parameter e, is an arbitrary function of this parameter only. The restriction 
of the group admitted by RG-manifold ( 0 ) on the solution of the BVP u = u(t,x,£), n = n(t,x,£ ) leads to zero 
equalities for two coordinates of the operator ( 0 ) in the canonical form - the conditions of functional self-similarity: 

r ? 1 + ^£uu x - i 2 u x - £ 3 u e = 0, rj 2 + ^e{nu) x - £, 2 n x - £ 3 n e = 0. (14) 

These equalities should be valid for any values of t, and certainly for t = 0, when dependencies u and n upon x are 
given by boundary conditions ©• This yields two linear relations between J 1 and J 4 : 

J 5 = U X J 1 , J 6 =N x J 1 + NJ 3 -N{U x ) u J 1 -£U x J 2 -NJ 3 ~NJ a . (15) 

Here, and in what follows functions U and N and their derivatives with respect to x should be expressed either in 
terms of u or in terms of y. Substituting ( 0 ) in (|l3]) gives the desired RG-symmetries with the RG-operator 

4 

r =r. R-i ■ ( ie ) 

i=1 


Ri 


X\ T 




(U x ) u 



J 1 


+ I £tU x — 


N 


J 1 — U J 1 


D + U X J 1 {£td x + d u ) , 


R 2 =X 2 + —J 2 D , R k = X k + —J k D, k = 3,4. 
n n 

We see that RG-symmetries for ( |ll|) , (jl^) are presented as a combination of symmetries of infinite-dimensional algebra 
with the infinitesimal operator (|13|). Any of the four operators R k (and their linear combinations with coefficients 
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that are arbitrary functions of e) contains the BVP solution u = u(t , x, e) and n = n(t, x, e) in the invariant manifold 
and enables to obtain group transformation of both group variables and different functionals of the solution (for a 
method of calculating transformation of a functional see Ref. fl7||). 

Generally, the renormalization group using is capable of improving a perturbation theory solution. As an example, 
consider the perturbative solution of (ll|), (jl2j) for small value of et 1 

u = U{x) — ( et)UU x + O (e 2 t 2 ) , n = N(x) — (et) ( UN X + NU X ) + O (e 2 t 2 ) . (17) 

This approximate solution in the limit (et) —> 0 is invariant with respect to RG transformation defined by the operator 
i ?2 with arbitrary eJ 2 y^ 0. Assuming J 2 = 1/e, we obtain the explicit expression for RG-operator 

R = - [(1 + £tU x ) (d t + eud x ) - eU x nd n ] , (18) 

n 

and invariance conditions written in the form of two first order DEs: 


ut + euu x = 0, (1 + etU x )(nt + eun x ) + enU x = 0. 


(19) 


Solving Lie equations which correspond to RG-operator 
enables to reconstruct the desired exact solution of ( 11 ), 


(|l8l ) (and coincide with characteristics equations for ©) 
( | 12 | ) from the perturbative solution © 


u = U(x — eut ), 


1 


n = 


1 + etU, 


-N(x — etu ), 


( 20 ) 


where U x should be expressed in terms of u. For example, in particular case of N(x) = N 0 exp(— x 2 ), U(x) = — x and 
e = 1/T the latter formulas describe the focusing of gaussian laser beam in geometrical optics 


T 

n = — - Nn exp 

T-t 



T 


t — T 


T 

u = x -, t < T. 

t-T ’ 


( 21 ) 


3.2. Now let us turn to a more general case of a 0. The Lie point symmetry group, admitted by RG-manifold 
is characterized by a canonical infinitesimal operator |l2| with six independent terms Xi = 1,..., 5 and X x 


5 

X = Xryc + ^ ( C , X , 
2=1 


/ + Ci fi ) d T + 





( 22 ) 


where coordinates fi and gi are linear combinations of r and x an d their first derivatives iq = (dr/dn) and Xi = 
(dx/dn) with coefficients depending only on v and n [^0|41|. For a particular case ip = 1 they are 


fi=T, 5i = x; f 2 = -(l/a)nxi, 92 = n; 

/ 3 = —r/2 + nn + (l/2a)nvxi, 53 = ~(v/2)t 1 + nxi] 

/ 4 = —(l/2)nx + vnn + [(l/4a)v 2 - n] nx i, (23) 

54 = (a/2 )t + (l/2)vx + vnx l + [an - (l/4)v 2 ] n; 

/.5 = (nn - r) - ar a , g 5 = nxl - ax a ■ 


” Evident” symmetries /i, 51 and / 2 , 52 describe dilations of r and x and translations along u-axis respectively for 
an arbitrary nonlinearity <p(n). Two more symmetries / 3 , 5 3 and / 4 , 54 appear due to a special form of the function 
ip = 1 under consideration. The symmetry /s, <75 involves the parameter a transformation along with transformations 
of dynamic variables. 

The operator X^ with coordinates / = ^(iqn), 5 = ^ 2 (v,n) that are arbitrary solutions of partial DEs 

£-(n/o)£= 0, £+£ = 0, (24) 

results from the linearity of basic Eqs. ©; it is an ideal of an infinite-dimensional Lie algebra formed by operators 
X \,..., A 5 and A/o* 
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The restriction of the group ( f23|) on the BVP solution means that coordinates / and g of the canonical operator 
(f 22 |) turns to zero on this solution, that is 


5 5 

f = ~E c ^ 9=-^Z c i9i- (25) 

2=1 2=1 

These relations express functions / and g in terms of /*, i = 1,..., 5 taken on a solution r = f(u, n), x — x( v i n ) 
of a BVP (exact or approximate). Substitution of (|2^) in ( p2| ) gives five RG-operators 

5 

R = V Cj(a)Ri , (26) 

2=1 


each being determined by corresponding coordinates /), and by a pair of functions A ®, B l 

Rl — (t — A 1 )^ + (x — B 1 )t}y , — A 2 ^-^ > 

R 3 = (-(r/2) - A 3 )d r - -B 3 9 x - (v/2)d v - nd n , 

i? 4 = (-(n/2) X - A 4 ) 5 r + ((a/2)r + («/2) X - S 4 ) (27) 

+ (—(l/4)n 2 + an) d v + vnd n , 

R 5 = ( t - A 5 ) <9r - i? 5 9x - nd n + ad a 

Here, ten functions A 1 , B’ are defined by expressions (H|) for /* and g 1 where one should replace r, % by f(n,v), 
X(n,v). Explicit formulas for RG-operators depend upon the specific solution of the BVP. For example, for the 
particular solution of the BVP (|)) with V = 0 and N(x) = cosh~ 2 (:r), described by @ 


r = 


( u /2) 1 / 2 


, 3/4 


(^ 


K 2 + 1 — K 


1/2 


1_n “4 ah 


1 (v/2\fa) x l 2 + (Vn 2 + 1 - k ) 

X = -oln 


1/2 


2 -(w/2i/a) 1 / 2 -I- (Vk 2 + 1 - k) 

functions A 5 , R 5 in (^) are expressed as follows fi(| : 

(c/2) 1/2 


1/2 


A 5 = — 


(v n 2 +1 — k) ^ + yrr^- 2 ^) 


R 3 = 


4a 3 / 4 -/1 + n 2 

{v/2) 1 ^ 2 (vAc 2 + 1 — k) ^ 

4a 1//4 -\/l + k 2 (Vk 2 + 1 — k — (v/2y/a)) 


STT^-3 K + 2^- V 

da 


(28) 


It should be noticed, that the solution of the presented above BVP is unique, but a number of RG-operators that 
give rise to this solution is different from one (in the first example case we have four RG-operators with arbitrary 
functions of (n, %), and in the second example five RG-operators with arbitrary functions of a). In the next section 
we will show that the number of RG operators may be enlarged to an arbitrary value, provided not only point but 
Lie-Backlund groups are taken into account. 


IV. RG AS LIE-BACKLUND SUBGROUP 

The method of constructing RG-symmetries from Lie point symmetries admitted by the original DE is naturally 
generalized to include Lie-Backlund (L-B) symmetries. The extension of the space of differential variables increases 
the amount of BVPs that allow restriction of a group on their solution. A complete set of RG-symmetries is obtained 
by appending L-B RG-symmetries to point RG-symmetries. In this section we present an example of constructing 
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L-B RG-symmetries of the second order for the BVP (®j). As in the previous section we use a transformed form of 
the basic equations ©• „ 

L-B symmetries admitted by the RG-manifold (|n^) are characterized by the same canonical infinitesimal operator 
(^||) where additional terms proportional to higher-order derivatives of r and x should be added in coordinates / 
and g. Similarly to first-order symmetries, these terms are linear combinations of r and \ and their derivatives 
n = (d l T/dn l ) and Xi = (9 l x/9n l ) with coefficients that depend only on v and n 39 -||]. For the second-order 
Lie-Backlund symmetries in a particular case = 1, we have five additional operators X- L with i = 7,..., 11 (the 
term with i = 6 corresponds to X,*, and is omitted in the sum, i.e. eg = 0 ) 


li 


li 


li 


X = ^ dXi = / + afi ) d T + g + ^ °i9i 9 X . 


(29) 


2=1 


2=1 


2=1 


It should be noted that expressions for all coordinates in 
recursive operators Ell Li , i = 1 , 2 , 3 


can be obtained by the action of the following three 


Li — 


Lr — 


D 


) — (n/ a)D n \ 

( 2nD n — 1 ( n/a)vD , 

, l 2 

= 

'n 0 j 

\ ~vD n 2 nD n 

2 nvD n 

n(v 2 /2a — 2n)D n — n \ 

-v 2 /2 + 2an)D n + a 

2 nvD n + v I 


(30) 


on the ’’trivial” operator with / = r and g = X (here, D n is the operator of total differentiation with respect to n). 
Below, we present only three of these five second-order L-B operators 


g$ = (l/2a)vxi + nr 2 + —nvx 2 ; 

Za 


(31) 


h = nr 2i g7 = xi + nx2\ 

fs = (l/2a)n [-Xi + vt 2 - 2 nx 2 \, 

/ 9 = (l/4)r - nn - (5/4a)vnxi + (~n + (l/4a)u 2 ) nr 2 - (l/a)vn 2 X 2 , 

gg = (3/4)uri — (2 n — (l/4a)u 2 ) xi + vnr 2 + (—n + (l/4a)i> 2 ) 2 • 

The procedure of restriction of the L-B group obtained on the solution of the BVP leads to expressions for f and q 
akin to © 


11 


11 


/=-E c ^ 9 = ~J2 


c i9i • 


(32) 


Substitution of (|32|) in ( p29| ) yields additional terms in the expression (26) for the RG-operator R that depends on 
higher-order derivatives of r and x 


r = J2 = E c ^ a ) ((/< - + (9i - s i )d x ). 


(33) 


i=1 


i=1 


Here functions A 1 and B l are given by the corresponding formulas for coordinates fi and gi to be evaluated on the 
solution f(n,v) and x(n, v). It appears that coordinates of L-B RG-operators are obtained from point RG-operators 
with the help of the above-mentioned recursive operators, hence, one can obtain L-B RG-operators of an arbitrary high 
order. Despite an unusual form, we still call them RG-operators since they possess the main property of RG-operators, 
namely, they contain a solution of the BVP in their invariant manifold. 

The procedure of using L-B RG-operators is not as simple as for point RG-operators. Yet we can describe two 
possible ways. 

Firstly, coordinates of canonical L-B RG-operators can be used to construct a set of relations, differential constraints, 
that are compatible with the original DEs and satisfy specific boundary conditions. The use of such constraints is 
described in section 6 . In the general case, for an arbitrary L-B group of a given order, coordinates of the corresponding 
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canonical operator can be treated as a set of differential expressions, zero equalities for which impose appropriate 
restrictions on the basic DEs, consistent either with physical or with symmetry conditions. These equalities can also 
be treated as embedding equations (see ^7j). 

Secondly, L-B RG-operators can be used to construct invariant solutions that automatically fit boundary conditions. 
It should be noticed that in some particular cases, L-B RG-svmmetries can be constructed from a L-B group with a 
finite number of operators. For example, RG-symmetry for (J9j) with boundary conditions V = 0 and N = cosh~ 2 (ir) 
appears as a linear combination of three L-B symmetries 

R = {fs + 2(iV + h))d T + (33 + 2(37 + g§))d x . (34) 

The desired solution of the BVP can be found as the invariant solution with respect to RG-operator ([dJ) and is 
presented by formulas (p8|). 

The recipe of constructing the L-B renormgroup formulated in this section goes far beyond a simple illustrative 
example for the BVP (||). In a similar way, L-B RG-operators are constructed for different BVPs of mathematical 
physics that admit L-B symmetries; other examples are presented in [^?| for the linear parabolic and modified Burgers 
equation. It is essential that when parameters entering into the equation and boundary conditions are involved in 
group transformations, coordinates of canonical L-B RG-operators contain not only first but higher-order derivatives 
with respect to these parameters. This means that in addition to recursive operators containing operators of total 
differentiation with respect to n (for BVP (§])), new recursive operators comprise operators of differentiation with 
respect to parameters, as well (oc D e and D a in the case of BVPs ([II]) and ©)• 

V. RG DEVISING BASED ON EMBEDDING EQUATIONS 

In this section we present a specific method of constructing RG-symmetries . which is based on embedding 

equations J3(|. It is of prime interest for physical systems described by ordinary differential equations (ODEs). In 
the context of the discussed model of quasi-Chaplygin media such equations arise, e.g., when constructing invariant 
solutions with respect to symmetries obtained. We demonstrate the idea of this method for the very simple BVP 

ut = f(t,u,a ); f = T, u = x. (35) 

Extension of the original differential manifold by adding, to the original equation, the embedding equation that 
appears as a linear first-order partial DE 


u r + f(T,x,a)u x = 0, (36) 

gives the desired RG-manifold, where u is now treated as the function of four variables {t, r, x , a}. Performing the 
group analysis for this RG-manifold involves boundary data and parameter a in group transformations, while the 
subsequent restriction of the group obtained on any solution of the BVP yields the desired RG-symmetries. We give 
two examples of such calculations for / = au 2 and / = u 2 + au 3 . 

5.1. In the event of / = au 2 the RG-manifold (p5|)-(|36l) is given by two equations 

u,t = au 2 , u T + ax 2 u x = 0 (37) 

that admit an infinite-dimensional Lie point algebra with five independent elements 

5 

X ^ ( Ot-iXi , 

2 = 1 

o o ( 38 ) 

X\ = d t + au 2 d u , X 2 = d T + ax 2 d x , 

X 3 = u 2 d u , X 4 = x 2 d x , X 5 = x 2 rd x + u 2 td u + d a . 

Here, functions aq and 012 depend upon five variables {t, r, x, a, zt}, whereas a*, i = 3,4,5 are arbitrary functions of 
three combinationsaf + ( 1 /zt), ar + {1/x), a. 

The procedure of restriction of the group obtained leads to the invariance condition 

U 2 (a 3 + aai + a 5 t ) — a\U t - a 2 U T — x 2 (a 4 + aa 2 + a 5 r)U x - a 5 U a = 0 (39) 
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to be fulfilled on an exact or approximate solution u = U(t, x , t, a) of the BVP 
the perturbative solution as an expansion in powers of a 

u = U(t, x, t, a) = x + ax 2 (t — r) + 0(a 2 ) , a<l. 


3); for example, one can take 


(40) 


Substituting ( }40| ) into (39) shows that the invariance condition ( |39| ) is fulfilled for 03 = 04 = a and arbitrary 01 , 0:2 
and 05 . Assuming a± = a.i = a = 0 and 05 = 1 in ( j38[ ) yields one of the RG-operators 

R = x 2 rd x + d a + u 2 td u , (41) 

which enables us to transform the perturbative solution of (EfHI) for small a <C 1 to the following exact solution 


1 — ax(t — t) 


This result is found by solving the Lie equations, that correspond to the RG-operator ([Il|). 

5.2. For another value of the function f = u 2 + au 3 , the RG-operator that is similar to © is given as follows 


R = (x 2 (l + ax)r + x) d x + (u 2 ( 1 + au)t + u) d u — ad a . 

The invarince condition for the solution of the BVP with respect to the RG-operator 
partial DE 


(42) 

has the form of the first-order 


— (ir 2 (l + ax)r + x) u x + au a + u 2 (l + au)t + u = 0 . 

Solving the characteristic equations for (^) (Lie equations) yields the following exact solution of the BVP 
f = u 2 + au 3 


(43) 

with 


1 1 

t — t — -ham 

x u 


x (1 
u (1 


au) 


+ ax) 


What all renormgroups obtained for the BVPs for the first-order ODE in the above examples have in common is 
that their operators depend upon arbitrary functions a;, which means that RG can be expressed in terms of different 
RG-operators with various particular expressions for their coordinates. This situation is the same as that one obtains 
for the BVP in the case of partial DE: different RG-operators yield the same unique specific solution of the given 
BVP contained in the invariant manifold of RG-operators. The previous procedure of RG constructing for the BVP 
for the ODE was based on the use of point groups. However, L-B groups can also be employed for constructing 
RG-symmetries for the first-order ODE, especially, in view of embedding equations (see Refs, in p7|). 

The structure of embedding equations depends not only on the form of the original equation, but also on the 
boundary conditions. This means that for given basic equations we may obtain different embedding equations. For 
example, if the function / in the r.h.s. of (35) depends upon x, 


Ut = f(t , x,a,u)\ t = r , u — x 

we arrive at the embedding equation 

Ur + f(r, x, a, x)u x = /(r, x, a, x) 


(44) 


dt 1 f x (t', x, a, u(t')) exp 


Z 

- J dt"f u (t",x,a,u(t")) 


(45) 


Hence, the RG manifold in this case is defined by a system of integro-differential equations ([S]) and ( |45| ) and one should 
employ the modern group analysis techniques which give a possibility of analyzing such equations, as well J|9],[50j . 


VI. RG AND DIFFERENTIAL CONSTRAINT 

In the previous section RG-manifold was obtained by combining an original DE and an embedding equation. 
More generally instead of an embedding equation, an additional differential constraint can be used that satisfy two 
conditions: firstly, it must be compatible with the original DE and, secondly, it should explicitly take boundary 
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conditions into account. This constraint naturally emerges when a coordinate of a canonical operator of the L-B 
RG admitted by BVP is assumed to be equal to zero. Adding this constraint to original equations we obtain the 
RG-manifold. 

6.1. To illustrate, consider first a BVP (|)j) with a = 0 which we rewrite using hodograph transformations in a 
simple form (compare with ©) 


Xn = 0 , Xv + T n = 0 . 


(46) 


L-B symmetries of this system of DEs are given by a canonical operator 


V = fd T + gd x , 


(47) 


with coordinates / and g depending upon v and derivatives r s + nxs+u Xs of an arbitrary order s > 0 


/ = F(v,Xs,t s ) - n 


dv + $> + i3 fs + Xk+id Xk ) 


k=0 


G, g = G(v,x s ,f a ), 


(48) 


T S =T S + nx s +i , T s = ( d s r/dv s ) , x* = ( d s x/dv s ) . 

Consider a particular case of a BVP (|9|) with boundary conditions defined by V(x) = —ex and arbitrary N(x). In 
terms of the variables r and x, these conditions are described, for example, by a pair of differential constraints 


X.VV — b 5 1~VV NvvXv NyXvV — 0 • 


(49) 


Here the dependence of N upon x is given in terms of v with the use of the above boundary condition. 

It is easily checked by direct substituting into (fl8]) that left-hand sides of these equalities are the corresponding 
coordinates g and / of the second-order L-B symmetry operator (|47|). Adding differential constraints ( |49| ) to the 
original equation (|l6|), we obtain the desired RG-manifold 

Xn — 0 j Xv ~b x n — 0 , Xvv — 0 > x vv N vv Xv — 0 • (^0) 


The latter admits a 17-parameter group of point transformations given by the following operators 

m 

X = Y^ CjXj , m = 17, 

i= 1 


(51) 


Xi = v 2 d v + v(2(n - N) + vN v )d n + (x{N - n) + rv)d T + vxd x , 

X 2 = vxd v + {x(n - N)+ v{xN v - r))d n + 2 rxd T + X 2 d\ , 

X 3 = -vdy +(N -n- v N v )d n , V 4 = vxd n - X 2 d T , 

X 5 = vd n , X 6 = (N - n) d n + x<9 x , X 7 = (n — N)d n + rd T , 

X 3 = d T , X g = vd T , Xio = (N — n)d T + vd x , 

Xn = d x , X r2 = xdr , Via = -v 2 d n + v\d T . 

V 14 = -d v - N v d n , Vi 5 = d n , Vi 6 = xd n , X 17 = -xd v + (r - x N v)d n ■ 

The usual procedure of restriction of the group obtained on a solution of the BVP ([it]) relates different coefficients in 
the sum ( |5l| ) and gives the desired RG operators 

13 

R=J2 c i R i . ( 52 ) 
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Ri = X \, i?2 = X2 , i ?3 = X3 + eXn , R4 = X4 , 
i? 5 = A5 + eXiQ , R e = Xq + eXyj , R-j = A7 + eX 16 , 

i?8 = Xg + £^ 15 ; R9 = Ag _ £ 2 -Yi6, i?io = -Y10 — £ 2 Ai7, 

-Rll = Xu + eXi4 , i?!2 = X12 + sX\Q , i?!3 = A13 . 

The exact solution of the BVP \ — ~ v / £ i T = (l/e)(n — N) is found either by solving Lie equations corresponding to 
any of these RG-operators, or as the intersection of all invariant manifolds. 

6.2. Now let us turn to a general case of a BVP (]|) with a/0. We shall consider the problem of constructing RG- 
symmetries using the RG-manifold given by basic equations in the form (nO f) an d the most simple differential constraint 
yielded by the linear combination of the second order L-B symmetry ([U]) /y, <77 and trivial infinite-dimensional 
symmetry f oc = 0, g 00 = -l 


ar v - nxn = 0 , Xv + T n = 0 , nr nn = 0 , Xn + n,\rm -1 = 0. (53) 

This differential constraint describes, in particular, a linear dependence of N upon x and V(x) = 0. The Lie point 
group admitted by the RG-manifold (^) is characterized by seven infinitesimal operators (use the formula (^TJ) for 
to = 7) 


X\ = —2vd v — 4 nd n — 6n(v/a)d T + (2% — 6 n + 3v 2 /a)d x , 

X 2 = vd v + 2 nd n + (t + 2 nv/a)d T + (2 n — v 2 /a)d x , X 3 = d v , 

X 4 = nd T — vd x , X 5 = (v/a)d T + In nd x , A" 6 = d T , A '7 = d x . 

The restriction of this group on the solution of the BVP with the above-mentioned boundary conditions leads to the 
three-parameter RG 


Ri — X\, i ?2 — X2, R3 — (1X3 + A4 . 

As in the previous case, the exact solution of the BVP r = nw, x = n — aw 2 /2 appears as an intersection of all 
invariant manifolds that correspond to these RG-operators. 

The characteristic feature of the described approach is the formulation of boundary data in the form of a differential 
constraint and the subsequent search of the group admitted by this constraint and basic equations. It is evident that 
there exists an infinite number of other differential constraints that adequately describe the same boundary data and 
the use of which leads to different RG algebras. As an example, we can point to differential constraints that arise 
from the zero equality of appropriate coordinates of the infinite L-B algebra. 

The example of RG-symmetry constructing on the basis of L-B symmetry reveals the practical importance of the 
latter and, on the other hand, demonstrates point symmetries that are not admitted by the original equation. The 
procedure of construction of RG-symmetries with the help of a differential constraint was also carried out in |37|] for 
the linear parabolic equation. 


VII. RG AS A SUBGROUP OF AN APPROXIMATE SYMMETRY GROUP 

An attractive method of RG constructing is that based on approximate symmetries pl[ . This method can be 
applied to systems described in terms of models based on DEs with small parameters. These small parameters allows 
us to consider a simple subsystem of the original DEs that usually admits an extended symmetry group inherited by 
the original DEs. Restricting this approximate group on the solution of the BVP yields the desired RG-symmetries. 
The merits of the described method is illustrated below for the BVP (^) with a small nonlinearitvparameter a 1. 
In terms of the variable w = v/a the following basic system of linear DEs is obtained instead of (|l0|) : 

r w - {n/tf(n))xn = 0 , Xw + ar n = 0 . (54) 

For a = 0, it admits an infinite L-B symmetry group 

X = fd T + gd x , (55) 
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characterized by an arbitrary dependence of the zero-order coordinates f = f° and g = g° upon n, r, x and the 
derivatives f s , \s of an arbitrary order 


f° = F° + J dw {(n/tp)Yg 0 } , g° = G° . 


(56) 


Here and below in (E 


Y = dn+ ^2 ( t s+i9t s + Xs+id Xs ) 


s=0 


d s r 

dn s 


Xs = 


d°X 

dn s ’ 


T s = T s W 


y- ( d p (n/V) 

V pJ dnP 


p=o 


Xs—p+l ) 


i 7 ' l (n,x s ,f s ) and G I (n,x s ,f s ) are arbitrary functions of their arguments, and expressions in curly brackets before 
integrating over w should be given in terms of t s . Xs, n, w. 

For 0 < a < 1, this symmetry is inherited as an approximate one by equations (|HJ) which thus represent an 
approximate RG-manifold. For example, for tp(n) = 1 the following result is obtained: 


f i = F* + Jdw | Zf- 1 + 'lYg i J , g i = G' + J dw {Zg l ~ l - Y f l ~ 1 } 

OO 

z = ^2 Ts + ld Xs > T s = T s - w{nXs+i + sXs ), i > 1 ■ 


(57) 


s=0 


One can see, that the symmetry of equations (|j) for a = 0 is inherited by the symmetry of these equations for a/0 
up to an arbitrary order of this parameter. It should be noticed that both zero-order and higher-order approximate 
symmetries may appear as Lie point symmetries or L-B symmetries, and this parameter may be involved in group 
transformations, as well. 

The restriction of the approximate group obtained on a particular solution of the BVP defines the specific form of 
the zero-order symmetries. It means that while constructing RG-symmetries for the BVP (^dj) in view of the bounary 
data from (j|), coordinates f°, g° and ’’integration constants” F l , G l , i > 1 are not arbitrary functions, but should 
be chosen so that relations f = 0, g = 0 satisfy desired boundary conditions r s =0, y = H{n) at w = 0. Provided 
that the functions F l and G l , i> 1 are also equal to zero in this case, boundary conditions are correlated with the 
form of functions /° and g°. In general, invariance conditions / = 0 and <7 = 0 appear as differential constraints (or 
algebraic relations) to be fitted by boundary data. 

Of special interest are such zero-order functions f° and g° for which infinite series ( |57| ) are truncated for some finite 
value of i = i m ax , and we arrive at finite sums. In this case, instead of an approximate group with respect to a small 
parameter a we obtain the exact symmetry group (compare with p7| , §11]). A simple example of this is given by the 
RG-operator (]34|). It is easily checked that in terms of n and w, the combinations of coordinates fz + 2(/y + fg) and 
g 3 + 2(gT + gg) are expressed as binomial in a, i.e. expressions for / and g are represented as zero-order and first-order 
terms f = f° + af 1 and g = g° + ag 1 , where /°, g° and f 1 , g 1 according to (p6|) and (57) are defined by the formulas 


/° = 2n(l - n)r 2 - nri - 2nw(xi + n\ 2 ), 
g° = 2n(l - n)x 2 + (2 - 3n)yi , 

/ 1 = t -nw 2 T 2 , g 1 = 2nwT 2 + wti + i( nw 2 Xi + ^Xi) ■ 


(58) 


^From here, in view of (^7|), it follows that higher-order corrections vanish, and we obtain an exact second-order L-B 
symmetry of DEs ( |54| ) at ip = 1 for arbitrary a /= 0; this symmetry gives rise to the exact solution ||38|,[39]] satisfying 
the boundary condition N = cosIR 2 (cc) defined by the zero order term g°. 

The arbitrariness in functions / , g° enables us to construct RG-symmetries for any boundary conditions. As an 
illustration, we present RG-symmetries for the BVP with 


H(n) = (ln(Vn )) 1 / 2 , 


(59) 
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describing space evolution (self-focusing) of the gaussian beam with the originally plane phase front at r = 0. To 
satisfy the initial distribution (|5f]), one can choose the following functions /° = 1 + 2n\Xi and 9 0 = 0. For this value 
of f° the inherited point group of the BVP is constructed with the help of formulas (|57|) and is given by the operator 

R = -2\d w + 2 ard n + ^1 + d T . (60) 

The invariance condition for the solution of the BVP with respect to RG with this operator is presented in the form 
of two partial DEs 

XXw ~ arxn = 0 , 2xt w - 2arr„ + 1 + (ar 2 /n) = 0 , 


the solution of which yields the desired approximate analytical solution of the problem 

2 


x 2 = {ant 2 — Inn) 1 — 2 Q{Vant' 2 ) 
Here the function Q(z) is expressed as follows 


x Q{V ant 2 ) 

v = -2 — - 
t 


1 — 2Q{V ant 2 ) 


(61) 


Q{z) = ze 2 ! 2 J dte * ^ 2 . 

o 


The first-order approximate symmetry obtained can be used to calculate a higher-order approximation of the RG- 
operator (|60|) and, thus, to improve the analytical solutio n (| 6 l| ). One can also obtain new-type RG-operators just by 
substituting the approximate solution ([n]) into formulas (|27[). 


VIII. CONCLUSION 


This paper presents a new approach to constructing RG-symmetries based on the mathematical apparatus of 
classical and modern group analysis. It differs from the traditionally used methods of constructing RGs in theoretical 
physics and is formulated as a sequence of the following steps: 

I) constructing the RG-manifold , that takes into account both basic equations and the corresponding boundary 
conditions; 

II) calculating the symmetry group , admitted by RG-manifold; 

III) restricting the group obtained on the solution of the B VP; 

IV) utilizing of RG-operators to find analytical expressions for solutions. 

As it was shown there exists a set of different algorithms for finding RG-symmetries. The choice of a particular one 
for a given physical problem depends on a mathematical model used for the problem description. 

It should be noted, that different methods of constructing RG-symmetries described above do not exhaust the 
suggested approach (see, e.g., |T2| . [~iTi|| ). Procedure of constructing RG-symmetries may combine different algorithms; 
for example, of interest is a simultaneous use of the method based on approximate symmetries and the invariant 
embedding method, and so on. 

Our approach reveals a close relation of functional self-similarity property (i.e., ” classical” RG-symmetry as an exact 
property of a solution) to an invariance condition of a BVP solution with respect to RG-operator. Mathematically, 
the latter is formulated as the vanishing condition for the coordinate of a canonical RG-operator on a solution of 


BVP. 

One can readily see that RG-operators may appear in the form, that is different from QFT case Jg), e.g., opera¬ 
tors of Lie-Backlund RG-symmetries. However, in some cases ”our” RG-operators can look like that ones in QFT 
renormalization group. For example, linear combination of operators a±X\ and 03(^3 + X 4 ) for the BVP (35) with 
oq = 1 , «3 = —a gives 


R = d t — ax 2 d x , 


(62) 


which is formally equivalent (with appropriate change of variables t = \nx, x = g and (3(g) = ag 2 ) to the differential 
operator for one-coupling massless QFT model in one-loop approximation. 

Up to now this approach is feasible for systems that can be described by DEs and is based on the formalism of 
modern group analysis. 
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It seems also possible to extend our approach on physical systems that are not described just by differential equa¬ 
tions. A chance of such extension is based on recent advances in group analysis of systems of integro-differential 
equations f|49 50 that allow transformations of both dynamical variables and functionals of a solution to be for¬ 
mulated H47| . More intriguing is the issue of a possibility of constructing a regular approach for more complicated 
systems, in particular to that ones having an infinite number of degrees of freedom. The formers can be represented 
in a compact form by functional integrals (or path integrals). 

This work was supported by Russian Foundation for Fundamental Research (project No 96-01-00195 and partially 
projects Nos 96-01-01297 and 96-15-96030). 
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